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1. INTRODUCTION 

The folding model has been used for years to calculate the nucleon-nucleus optical 
potential (see, e.g., Refs. |]-[5]|) and inelastic form factors (see, e.g., Refs. [0,0]). It can 
be seen from the basic folding formulas that this model generates the first-order term of 
the microscopic optical potential that is derived from Feshbach's theory of nuclear reac- 
tions ||. The success of this approach in describing the observed nucleon-nucleus elastic 
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scattering data for many targets suggests that the first-order term of the microscopic op- 
tical potential is indeed the dominant part of the nucleon optical potential. In the same 
way, the inelastic (folded) form factor is also the most important input for the analysis 
of inelastic scattering data within the distorted wave Born approximation (DWBA) or 
coupled-channel approaches. 

The basic inputs for a single-folding calculation of the nucleon-nucleus potential are 
the nuclear densities of the target and the effective nucleon-nucleon (NN) interaction. 
If one has a well tested, realistic effective NN interaction, the folding model is a very 
useful approach to check the target nuclear densities. A popular choice for the effective 
NN interaction has been one of the M3Y interactions which were designed to reproduce 
the G-matrix elements of the Reid || and Paris |L0| NN potentials in an oscillator basis. 
Although these density -independent M3Y interactions were originally developed for use in 
the DWBA analysis of (p,p') reaction, they have been used much more often in the double- 



folding calculation of the heavy-ion interaction potential at low and medium energies [11 



Intensive studies of the refractive nucleus-nucleus scattering during the last decade has 
shown that the simple M3Y-type interaction failed to give a good description of the data, 
and this has motivated the inclusion into the original M3Y interactions of an explicit 
density dependence [I^-CT, to account for the reduction in the (attractive) strength of 



the effective NN interaction that occurs as the density of the nuclear medium increases 
[llj. With parameter values chosen to reproduce the observed nuclear matter saturation 
density and binding energy within the Hartree-Fock (HF) approach, the new density 
-dependent M3Y interactions have been carefully tested in the folding analysis of the 
refractive nucleus-nucleus elastic scattering |l2|-p^]. In the same HF scheme, parameters 
of the density dependence of the M3Y interaction are directly associated with the nuclear 
incompressibility K, and it has been shown that K values ranging from 240 to 270 MeV 
are the most appropriate for the cold nuclear matter |15| (these values are rather consistent 
with what can be extracted from the calculations of the isoscalar giant monopole resonance 
within the relativistic mean field framework; if non-relativistic Skyrme or Gogny forces 
are employed, slightly lower K values are obtained [|T^]). The isospin dependence of this 
interaction was also tested in a HF study of the asymmetric nuclear matter and in the 



calculation of the heavy- ion potential between nuclei with non-zero isospins [fjj|. It is 
therefore of interest to use the new interaction in the folding calculation of the proton- 
nucleus potential for the analysis of the elastic and inelastic scattering of exotic nuclei on 
proton target. 



In the present work, we adapt the most recent version of the double-folding model [19 



for the nucleus-nucleus potential to a single-folding formalism for the elastic and inelastic 
nucleon-nucleus potentials, using the density- and isospin dependent M3Y interaction. 
By using the nucleon-nucleus elastic and transition potentials folded with the nuclear 



densities from a self-consistent microscopic model [20], we have performed a detailed 
folding analysis of the recent elastic and inelastic 30,32 S+p scattering data measured at 
E lah = 53A MeV |21] in GANIL, as well as 34,36,38,40g +p data measured earlier 

at .Biab = 30, 28, 39 and 30 A MeV, respectively. The contributions from isoscalar and 
isovector parts of the proton-nucleus optical potential and inelastic form factors were 
treated explicitly in each case to study the isovector mixing effect in the proton-nucleus 
scattering as one goes along the isotopic chain, from the proton rich 30 S to the neutron 
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rich (short lived) 38 > 40 S isotopes. 

2. THE FOLDING MODEL 

2.1. General formalism 

In the folding model, the proton-nucleus potential V is evaluated as a Hartree-Fock-type 
potential 

Vi = Y,[< ijMij > + < ij\vEx\ji >], (!) 
jeA 

where Vp and vex are the direct and exchange parts of the effective NN interaction be- 
tween the incident proton % and nucleon j in the target A (in the case of A scattering on 
proton, A is still treated in our formalism as 'target' as given by the inverse kinematics). 
The antisymmetrization of the system is done by taking into account the knock-on ex- 
change effects (the interchange of nucleons i and j). Due to the antisymmetrization, the 
potential is, in general, nonlocal in coordinate space. An accurate local approximation is 
usually obtained by treating the relative motion locally as a plane wave [0,0, and one can 
reduce the energy-dependent (central) proton-nucleus potential ([!]) to the following local 
form 

V(E, R) = J [p(r)v D (p, E, s) + p(R, r)v EX (p, E, s)j (k(E, R)s)} d\ (2) 

where s = r — R, R is the vector joining the center-of-mass of the target and the incident 
proton, p(r, r') is the nonlocal (one-body) density matrix (DM) for the target nucleus 
with p(r) = p(r, r) and k(E, R) is the local momentum of relative motion determined as 

k\E,R) = %[E c . m . - V(E,R) - V C (R)}. (3) 
n 

Here, p is the nucleon reduced mass, V(E, R) and Vq{R) are, respectively, the total nuclear 
and Coulomb potentials evaluated in the entrance (elastic) channel. The calculation of 
the localized exchange part of the nucleon-nucleus potential [second term in Eq. (H)] still 
contains a self-consistency problem and involves an explicit integration over the nonlocal 
nuclear DM. Normally, the density p(r) is taken either from a nuclear structure model or 
directly from electron scattering data. Therefore, the calculation of the exchange potential 
in Eq. (|2]) is usually done by using a realistic approximation for the nonlocal DM p],|26| 

p(R,R+s) ~ p (R + l)j 1 (k F (\R+~\)s 

where ji(x) = 3ji(x)/x = 3(sinx — xcosx)/x 3 . (4) 

To accelerate the convergence of the DM expansion, Campi and Bouyssy j2f| have sug- 
gested to choose the local Fermi momentum fcj?(r) in the following form 



k F (r) 



r(r) - \v 2 p{r) 



1/2 

(5) 



I 3p(?" 

Assuming this prescription, we choose further the extended Thomas-Fermi approximation 



27| , |28| for the kinetic energy density r(r) in the evaluation of the local Fermi momentum 
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kp(r). Note that the local approximations made in our approach are essentially the same 
as those adopted earlier in the folding calculations of nucleon-nucleus optical potential 
and inelastic form factor ||[7| . 
In this connection, it should be noted that there exists a more sophisticated version of 
the single-folding model |5| where the nonlocal exchange potential is treated exactly in 
the Schrodinger equation for the scattering wave function. In this approach one calculates 
the nonlocal nucleon-nucleus potential using the explicit expression for each single-particle 
wave function \j > taken, e.g., from the shell model. Therefore, this rigorous approach 
cannot be used in a general case, when the target wave function is simply represented by 
a local density distribution p(r). 

2.2. Explicit treatment of the isospin dependence 

Exotic nuclei usually have non-zero isospin and it is necessary to make explicit the 
isospin degrees of freedom. The spin-isospin decomposition of the (central) NN interaction 
is 

vd(ex) (p, E, s) = Vo {EX) (p, E, s) + Vi {EX) (p, E, s) {tr.tr') 

+ v° {EX \ P ,E,s)(t.t') + v° {EX \p,E,s)(ct. ( t'Kt.t'). (6) 

The contribution from the spin dependent terms (v w and t>n) in Eq. (|) to the central 
nucleon-nucleus potential (|2p is exactly zero for a spin-saturated target. Even for an odd 
nucleus, this contribution is at most of A' 1 effect [^9| and is usually neglected in the 
folding calculation of the central nucleon-nucleus potential. 

By writing the nuclear densities in Eq. (|2|) explicitly in terms of the proton (p p ) and 
neutron (p n ) densities, one can represent the proton-nucleus potential (|) in terms of 
isoscalar (V IS ) and isovector (V IV ) parts 

V(E, R) = V IS (E, R) + V IV (E, R). (7) 

Each term in Eq. ([?[) has contributions from both the direct and exchange potentials 

V IS (E,R) = j{\p p {r)+ Pn {r)]v°{p,E,s) 

+ [p p (R, r) + p n (R, r)]vg x {p, E, s) Jo (k(E, R)s)}d 3 r. (8) 



V IV (E,R) = J{[p p (r)-p n (r))v°(p,E,s) 

+ \p p (R, r) - p n (R, r)]^ x (p, E, s)j (k(E, R)s)}d 3 r. (9) 

One can see that the V IV term (microscopic form of the symmetry or Lane potential) is 
entirely determined by the difference between the proton and neutron densities. 

By using the local approximation (|j) for the nonlocal proton and neutron (r = p, n) 
density matrices 

p T (R,R + s) ~ ^(fl+^/^KIH+fl)*) = /rOR + f), 

h .r M fro 2 f , 5C s [Vp T (r)] 2 5VV(r) ) 1/2 

where k F {r) = \ 3vr p r r + + \ , (10) 

[ J 3p^(r) 36p r (r) J 
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the proton-nucleus potential ©— © can be obtained as 

V IS{IV \E,R) = J{p p (R+s)±p n (R+s)}v° m (p,E,s)d 3 s 

) (p,E,s) 3o (k(E,R)s)d 3 s. (11) 



+ 



f p (R+-)±f n (R 



2' 



W M(01) 1 



In Eq. (0), Cs is the strength of the so-called Weizsacker term representing the surface 
contribution to the kinetic energy density. We have adopted the commonly accepted 
value Cs^l/36 |27J which ensures a fast convergence of the density matrix expansion. 
The local Fermi momentum k T F (r) is evaluated using the ground state densities only. The 
local approximation ( [TDD was shown in a recent folding analysis to be at the level of 
1% accuracy. 

2.3. Effective density dependent NN interaction 

The recently parameterized CDM3Y6 interaction [|lj| , based on the G-matrix elements 
of the Paris NN potential ||10|| , is used in the present folding calculation. The energy- and 
density dependences are factorized out as 

v oom( E >Pi s ) = 9(E)F(p)v^\s). 
The explicit density dependence was introduced in Ref. [|1^ 



(12) 



F(p) =C[l+aexp(-/?p)- 7 p], (13) 
with parameters adjusted to reproduce saturation properties of nuclear matter and to 



yield a nuclear incompressibility K = 252 MeV in the HF approximation fl2| , |T5 
C = 0.2658, a = 3.8033, (3 = 1.4099 fm 3 , 7 = 4.0 fm 3 . 



(14) 



The 'intrinsic' energy dependence of the interaction (to be expected if one regards the 
effective NN interaction as representing a reaction- or G-matrix of the Brueckner type 
[]3 Of) is contained in the factor g(E) « 1 - 0.0026e, where e is the bombarding energy per 
nucleon (in MeV). The radial strengths of the isoscalar and isovector components of the 
central M3Y-Paris interaction [IT0I1 can be obtained ITS! in terms of three Yukawas 



D(EX) 
00(01) ^ 



d{ex) ( .exp(-R u s) 



00(01) 



R u s 



(15) 



where the explicit Yukawa strengths are tabulated in Table [I]. 

In nucleon-nucleus scattering the most important interaction induced by the nucleon 
spin is the spin-orbit coupling which is present in both elastic and inelastic channels. The 
spin-orbit potential arises naturally in the folding model if the effective NN interaction 
itself has a two-body spin-orbit term 



v LS (s)L.S = v LS (s) -[(r< 



(16) 



For simplicity, we assume that the spin-orbit part of the CDM3Y6 interaction has the 
same density- and energy dependences as the central part ( |12"D 



v LS (E,p, s) = g{E)F{p)v LS {s) 



(17) 
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Table 1 

Yukawa strengths of the central 
tion pl|§- 



lq) and spin-orbit (Qq) components of the M3Y-Paris interac- 



V 


R u 
(fm- 1 ) 


(MeV) 


(MeV) 


(MeV) 


Y£ X {v) 
(MeV) 


Yi%) 
(MeV) 


(MeV) 


1 


4.0 


11061.625 


313.625 


-1524.25 


-4118.0 


-5101.0 


-1897.0 


2 


2.5 


-2537.5 


223.5 


-518.75 


1054.75 


-337.0 


-632.0 


3 


0.7072 


0.0 


0.0 


-7.8474 


2.6157 


0.0 


0.0 



The radial strength of the spin-orbit components (with the total isospin T=0 and T=l) 
of the M3Y- Paris interaction [JToJ] can also be obtained in terms of Yukawas 



«£?(») = E^M^tr^. < 18 > 

with the explicit Yukawa strengths tabulated in Table [I]. 

2.4. Multipole decomposition 

For a consistent description of the elastic and inelastic scattering by the folding potential 
we need to take into account explicitly the multipole decomposition of the nuclear density 
nj,E01 that enters the folding calculation (^)-flTTD 



Pr(r) = E< JaMaMJa'Ma, > C x pl(r)[i x Y Xll (r)Y, where r = p,n; (19) 

X/j, 

J a and J a' are the target spins in the ground state and excited state, respectively. Usually, 
the A = term represents the ground state density (monopole excitation is a special case 
and not considered here) and a single multipole A ^ dominates in the transition to 
an excited state with spin J' A . In such a case, the corresponding term in the sum (|T9|) 
represents the nuclear transition density for the excited state. Following Satchler and 



Love [IlJ, we have chosen here a normalization such that C = VAn and C\=l for A ^ 0. 



The neutron and proton transition moments are further determined as 

POO 

Ml = / p T x (r)r x+2 dr. (20) 
Jo 

We adopt the same definition of the reduced matrix element as that by Brink and Satchler 
|31"H , and the nuclear transition density is such that the reduced electric transition rate for 
a 2 A -pole excitation is obtained from the proton transition moment as B(EX f ) = e 2 |M^| 2 . 

The corresponding multipole decomposition of the folded potential (|7|) can then be 
written as 

V{E,R) =£< JaMaMJa>Ma' > C x [v x IS (E,R) + V X IV (E, i?)J [i x Y^(R)]*. (21) 

Xfi 

The isoscalar and isovector parts of the central folded potential ( pip consist of the corre- 
sponding direct and exchange components 

Vl S{IV \E, R) = V^ S{IV \\, E, R) + Vi s x {iv \\, E, R). (22) 



7 



We note that if one writes the multipole expansion of the density ([19]) explicitly through 
the proton and neutron deformation parameters /3 T , one would end up with the multipole 
components of the folded potential that depend upon those (3 values. Such a version 
of the folding model will be suitable for the analysis of inelastic proton scattering from 
deformed nuclei, as has been studied earlier by Hamilton and Mackintosh in their density 
dependent folding model |32 . 



2.5. Direct potential 

Using the folding formulas in momentum space fll| the central direct potential can 
be calculated simply with the density dependent M3Y interaction (0). For the elastic 
scattering we have Ja' = J a an d A = 0, and the density components Pq{t) 111 ^q- (HD 
are just the proton and neutron ground state densities. Denoting the total ground state 
density as po(r) = Po( r ) + Po( r )i the direct part of the central elastic potential (|22|) can 
be obtained in the following form 



V^ IV \X = 0,E,R) = 9 -^± f A^ IV \ q )vg m ( q )j ( q R) q 2 dq, 



(23) 



where the Fourier transforms of the direct interaction and ground state density profile are 



(?) 



A, 



^00(01) 1 

IS(IV) 




47T 



47T 



Voo(oi)(r)j (qr)r 2 dr, 

[Po( r ) ± Po ( r )] F(p (r))j (qr)r 2 dr. 



(24) 



In evaluating the inelastic form factor (or transition potential) one needs to include the 
medium corrections implied by the use of a density dependent NN interaction (see, e.g., 
Refs. [fF| |I~9"| , |3~3f ) . In a consistent or 'dynamic ' treatment of the density dependence, the 
change in the density due to the excitation, p — > po + Ap, will also change the effective 



NN interaction 



F(p)v D{EX )(p, s) 



F(po) + Ap 



dF(p ) 



dpo 



vd(ex){s) 



(25) 



This prescription happens to be exact for the excitation of a single phonon of a 2 A -pole 

In this case, the direct transition potential fl2T 



harmonic shape vibration j34 
2 A -pole excitation of the target can be written in the following form 



for a 



V^ IV \X,E,R) 



9(E) 
2tt 2 



A{ S (/V) (g) + B^ lv \q) v»(q)j x (qR)q 2 dq, (26) 



} IS(IV). 



00(01)' 



where the Fourier transforms of the transition density profiles are 



.4 



and B 



IS {IV) 
A 

IS (IV) 
X 



(<?) 
(<?) 



47T 



Atx 



[p\{ r ) ± PaO)] F(p (r))j x (qr)r 2 dr, 
[Po» ± P n o(r)} Px(r) ^^- jx(qrydr. 



(27) 



Here, the total transition density is denoted as pxij') = PaM + Px( r )- ft i s eas y to see 



that5f (/F) (g) comes from the 'dynamic 7 treatment of the density dependent interaction 
discussed in Eq. (|25|). 



s 



2.6. Exchange potential 

The self-consistent (local) exchange potential has to be calculated by an iterative pro- 
cedure [|nj, and the exchange part of the elastic potential ( |22] ) can be evaluated as 



V^ IV) {\ = 0,E 7 R) = 2^(E)^Gf (/y H^^)<ofoi)(s)jo(M^^)^)^, (28) 
where Gf {IV) (R, s) = £ [f^y(x), s) ± ft(y(x), s)} F(p (y(x))dx, 



s 2 



y{x) = JR2 + - + R SX and ft(y, s) = pl(y) 3l (k F (y)s) . (29) 



After some integral transformation using the expansion formula |35| for spherical har- 
monics of the mixed argument Y\n(R + s), the exchange transition potential ( p2|) for a 
2 A -pole excitation of the target can be written as 

v£™(\,E,R) = 2ME) r[Gi s ^(R,s) + Hi s ^(R,s 

Jo 1 

x v^ 01) (s)j (k(E,R)s)s 2 ds, (30) 
where the exchange kernels are 



G^ IV \R,s) = fl* £ [fMx), s) ± K(y(x), s)] Fi ^ ]) dx, 
H^ IV \R,s) = R X jy p (y(x),s)±f-(y(x),s)]Pf^ 



dF(p (y(x))) 



dpo(y(x)) 



Gfo(3l) 



with fl{y, s) = p\(y)ji (kp(y)s) ■ Note that we have implied the same dynamic treatment 
of the density dependence ([25|) which leads to the hI S( * IV \r,s) kernel. 

2.7. Spin-orbit and Coulomb potentials 

In the elastic channel, in addition to the central potential V IS( * (E, R), the nucleon- 
nucleus optical potential has also the spin-orbit part V^A = 0, E, R)(l.a). In many 
cases, a phenomenological Thomas form for the spin-orbit potential suffices for a good 
description of the elastic data if the strength is adjusted properly. Within the folding 
model, the spin-orbit potential can be evaluated microscopically using the two-body spin- 
orbit NN interaction and the nuclear density of the target. We have chosen the local 
approximation developed by Brieva and Rook || and evaluated Vj,s(A = 0, E, R) using the 
spin-orbit component of the CDM3Y6 interaction flTTD-(|T%D, with the isospin dependence 
treated explicitly, 



pV ' J R dR v ' J R dR 



(32) 



%(E,R) = J Q C<} vi ) s (s)[l+j 1 (k(E,R)s)}s 4 ds, 

$ n (E,R) = ~ r{v^(s)[l+Uk(E,R)s)]+v^(s)[l-j\(k(E,R 
2 Jo 
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In its general formulation, the spin-orbit interaction also contributes to the inelastic scat- 
tering channel and one needs to explicitly calculate the scattering amplitude in the DWBA 
formalism, using the spherical tensor expansion of the spin-orbit interaction and explicit 
wave function for each single-particle configuration involved in the excitation |]36| , |38| . 
Such a method is not applicable if the target excitation is simply represented by a nu- 
clear transition density p\(r). Here, we have adopted the same local approximation || as 
in the elastic case and inserted ([19]) into the local density expression. In the first-order 
approximation, the spin-orbit transition potential can be obtained as 

where & T {E, R) is determined by the same formula fl3~3"|). Note that for the natural-parity 
excitations considered in this work, the central transition potential gives a dominant 
contribution to the inelastic cross section, and the transition spin-orbit potential ( |34] ) has 
only a minor effect. 

The proton-nucleus optical potential or inelastic transition potential, calculated by the 
folding model, has to be supplemented by a corresponding Coulomb potential. It is 
straightforward to use the same folding method to evaluate microscopically the proton- 
nucleus Coulomb potential, using the (target) charge density matrix 

2 

V C {E, R) — [ | 6 [pchargc(r) - Pchar ge (i*, r)j (k(E, R)\r - R\)\ d 3 r. (35) 
J \r — Jrt\ 

One can use the same local approximation (|n]) for the nonlocal proton density matrix and 
then take into account the finite proton size|TTJ to obtain the charge distribution PchargeM 
for the calculation of the Coulomb potential (|3~5"D . Vc(E, R) can then be expanded into a 
multipole series analogous to that of Eq. fl2~T|) for the nuclear potential. 

In the optical model (OM) analysis of elastic proton-nucleus scattering, the Coulomb 
potential Vc(X = 0, E, R) is usually represented by the Coulomb potential between a 
point charge and a uniform charge distribution of radius Rc = rcA 1 ^ 3 . This option 
for the elastic Coulomb potential can be shown to have about the same strength and 
shape at the surface as the microscopic potential obtained from Eq. (|35|). For convenience 
in referring or comparing with other results, we have chosen this option for the elastic 
Coulomb potential in our OM calculation. 

In the DWBA analysis of inelastic proton-nucleus scattering with a 2 A -pole excitation 
of the target, the transition Coulomb potential Vc(X, E, R) has been often taken in a 
model- and energy independent asymptotic form that can be expressed in terms of the 
reduced electric transition rate B(EX) (see, e.g., Refs. fl9| , |29"|j ). Since a correct input 
for the Coulomb form factor is substantial in the analysis of the low-lying electric type 
excitations, we have used in the present work the microscopic Coulomb form factor given 
by Eq. (0). 

2.8. Complex elastic and inelastic potentials 

The original M3Y interaction (|15"D is real, and the formalism presented above can be 
used to generate the real parts of the elastic and inelastic nucleon-nucleus potentials only. 
These must be supplemented by imaginary potentials which account for the absorption 



p(E,R) — 



dR 



R dR 



(34) 
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into other channels that are allowed energetically but are not considered explicitly. This 
absorption contributes imaginary parts to both the diagonal (elastic) and the off-diagonal 
(inelastic) potentials. It is common to resort to a hybrid approach by using the folding 
model to generate the real part but to use the standard Woods-Saxon (WS) potential for 
the imaginary part, and the local nucleon-nucleus optical potential is 

U (E, R) = V Q {E, R) + iW (E, R) + V LS (X = 0, E, R)(l.<r), (36) 

where the real and imaginary parts of the (central) optical potential are given by 

V {E,R) = V c (R) + N R [V D (0,E,R) + V EX (0,E,R)}, 

w r E = Wv 4iy 5 exp((r - R w )/a w ) 

1 + exp {{R — Rw)/a>w) [1 + exp((r - R w )/a w )] 2 

Here Vc(R) is the Coulomb potential between a point charge and a uniform charge distri- 
bution of radius Rc = rcA 1 / 3 , R\y = rwA 1 / 3 and the real optical potential Vd(ex)(0, E, R) 
is calculated by using Eqs. ( |23"1) and (^8|). The choice of the parameters for the WS 
imaginary potential can be made using the available systematics for the nucleon optical 



potential. The accurate CH89 global systematics p9| , based on several thousands data 
sets for the nucleon-nucleus elastic scattering at energies ranging from 10 to 65 MeV, 
is used in this work to parameterize the imaginary part of the optical potential. The 
normalization factor Nr for the real (folded) optical potential is adjusted in each case to 
fit the elastic data. In the present folding approach, Nr is an approximate way to make 
small adjustments that may be needed to take into account the higher order contributions 
to the (real) microscopic optical potential, the so-called 'dynamic polarization potential' 
(DPP) in the Feshbach's formalism ||. It is obvious that value of Nr should remain close 
to unity for this procedure to be reasonable. 

We stress that a reaction- or G-matrix calculation for a single nucleon incident on nu- 
clear matter []30| can lead to a complex effective NN interaction (like the JLM effective 
interaction), where the absorption is associated with the finite mean free path of nucle- 
ons in nuclear matter. The imaginary part of such an effective NN interaction (obtained 
in a standard local density approximation) does not account in principle for the impor- 
tant sources of absorption in finite nuclei like the excitation of surface modes, nucleon 
transfer and breakup reactions, etc. These non-elastic processes contribute to the energy 
dependent, nonlocal and complex DPP that gives rise to the imaginary part of the optical 
potential 0. 

We apply further the same hybrid scheme to the inelastic (transition) potential, and 
the transition potential for a 2 A -pole excitation is 

U X (E, R) = V X (E, R) + tW x (E, R) + V LS (X, E, R)(l.tr), (38) 

where the real and imaginary parts of the (central) inelastic potential are given by 

V X (E,R) = V C (X,E,R)+V D (X,E,R) + V EX (X,E,R), 

Wx(E , R) = (39) 

Thus, the real nuclear, Coulomb and spin-orbit transition form factors are calculated by 
the folding approach, while a conventional approach of the collective vibration model is 
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used to obtain the imaginary transition form factor. In the present work we have assumed 
5\ equal to the charge deformation length, which is determined from the reduced electric 
transition rate B(EX f) |29|] as 

4tiJB(EX T)/e 2 , 1 ! pl(r)r x+l dr 
S x = V where < r A-i > J ^40) 



The obtained charge deformation lengths for the lowest 2 + states in Sulfur isotopes are 
given in Table |2|. We note that the real central and spin-orbit inelastic form factors in 
Eq. ( j39|) were used in the DWBA calculation as given by the folding model, without a 
normalization factor. In the past, one has often used the same normalization factor Nr 
for both the elastic and inelastic folded (central) potentials. Such a method is 'consis- 
tent' if one considers Nr of the real (folded) optical potential as scaling factor of the 
effective NN interaction for a particular target at the given incident energy. In our case, 
parameters of the CDM3Y6 interaction have been fine-tuned, in the HF calculation, to 
reproduce the saturation properties of nuclear matter (with a realistic value of the nuclear 
incompressibility K « 252 MeV ||15|| ) and energy dependence of the nucleon- nucleus op- 



tical potential |T2| , |37| . Thus, the CDM3Y6 interaction should be suitable for any target 
(at different energies) and no further scaling of the interaction is expected (because it 
would affect the established HF results, on which the interaction is based). Therefore, 
one needs to consider the normalization factor Nr of the real (folded) optical potential 
in ([37] ) approximately as an effect caused by contribution of the DPP to the (real) micro- 
scopic optical potential only. Given the effective NN interaction carefully parameterized, 
it is reasonable to use the folded inelastic form factor without re-normalization in the 
(one-channel) DWBA calculation of inelastic scattering if the higher-order contribution 
from the DPP to the real inelastic folded form factor is weak. We will discuss this aspect 
again in Sec. | 

All the OM and DWBA analyses of elastic and inelastic proton scattering were made 
using the code ECIS97 written by J. Raynal [j3~8fl . 

3. THE NUCLEAR DENSITIES 

Different kinds of nuclear densities can be used as input of the present folding approach. 
Actually, Eqs. (|1|) and (pC|) refer to a macroscopic description (see, e.g., p. 579 of Ref. |29| ) . 
and the expression (^) should be interpreted as a generic deformed density of an excited 
spherical nucleus, in which the A 7^ multipole terms represent the difference with respect 
to the ground state density and may be used as transition densities for inelastic transi- 
tion to natural parity states of multipolarity A, provided they are properly normalized. 
On the other hand, the present folding approach can naturally accommodate, within its 
framework, the microscopic nuclear ground state and transition densities (see also p. 657 
of Ref. p9| ) and we believe it is worth to study the predictive power of the model by 
using these microscopic nuclear structure inputs. In this section, we briefly describe the 
structure model used to obtain the nuclear ground state and transition densities, having 
in mind the application of the model to the neutron rich Sulfur isotopes and giving some 
details relevant to this case. 
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We work in the framework of self-consistent calculations based on Hartree-Fock plus 
Random Phase Approximation (HF+RPA) with the use of Skyrme interactions. This 
nuclear structure model has been used for a long time to determine the structure prop- 
erties of the nuclear vibrational states, especially, the giant resonances. Recently it has 
been extended to include pairing correlations by performing quasi-particle RPA (QRPA) 



calculations on top of HF+BCS fl20| , |21f (see also Ref. ||40| ). In the study of nuclei along 
extended isotopic chains the inclusion of pairing is indeed called for, and its effect is 
expected to be important in the open shell isotopes, especially in the case of low-lying 
excitations. 

3.1. The nuclear ground state densities 

In this work, we perform spherical HF+BCS calculations in coordinate space using the 
SGII parameterization ||4j of the Skyrme interaction. We choose a simple prescription 
for the pairing interaction by assuming a constant pairing gap given by A = 12/a/A MeV 



i5fl . The use of a constant pairing gap produces unrealistic results (states at relatively 
high energy acquire quite large occupation probabilities) unless a cutoff is set in the single- 
particle space such that states above this cutoff do not feel any pairing interaction. In 
the case of the Sulfur isotopes, this cutoff is chosen such as to include all the subshells of 
the major shell to which the Fermi level belongs. The ground state neutron and proton 
densities are then obtained as 

P T o(r)=J2vlMr), (41) 

a 

where y9 Q (r) and v a are the radial wave function and BCS amplitude of the quasi-particle 
state \a >, respectively. The summation on a runs over neutrons (protons) for r=n 
(r=v)- 

3.2. The nuclear transition densities 

We briefly recall the main steps of the QRPA calculations p0 |. Based on the HF+BCS 
quasi-particle basis, the QRPA equations are solved in the configuration space, using 
the standard matrix form. The residual interaction between quasi-particles is derived 
from the Skyrme force used in the HF+BCS calculation, without including the pairing 
contribution. The continuous part of the single quasi-particle spectrum is discretized 
by diagonalizing the HF+BCS Hamiltonian on a harmonic oscillator basis. The size of 
the QRPA model space, i.e., the number of two quasi-particle configurations included, is 
chosen large enough to exhaust the appropriate energy- weighted sum rule. Solving the 
QRPA equations provides us with the energies uj k as well as the wave functions of the 
excited states \k > of a given multipolarity A. These wave functions are normally given 
in terms of the well-known amplitudes X£g and Y£ g , with af3 denoting the involved two 
quasi-particle configurations. The nuclear transition to the excited state \k > can be 
characterized by the corresponding (local) proton or neutron transition densities Ap T K (r) 
defined as 

Ap T K {r)=< K \'£S{r-r i )\g.B.>, (42) 

i 

where the summation on i runs over neutrons (protons) for r=n (r=p). 
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Table 2 

The excitation energy u, reduced electric transition probability B{E2 \) and ratio of quadrupole 
transition moments M = (M^+ / M^ + ) / (N / Z) for the lowest 2 + states in Sulfur isotopes, as given 
by QRPA (see Sec. |||). The 'experimental' M exp ratio is obtained by scaling the proton part 
of the QRPA transition density to reproduce the experimental B(E2) exp value and adjusting 
the neutron part to the best DWBA fit to the inelastic scattering data. The 'phenomenological' 
-^Phenom ratio was obtained in Refs. [^J2^] based on standard collective model form factor. The 
experimental values for uj exp and B(E2) exp are taken from Refs. [41-43]. The charge deformation 



lengths S2 are determined from the corresponding B(E2) cxp values using Eq. 



Nucleus 


WQRPA 




B (E2) qrpa 


B(E2) exp 


62 


Mqrpa 


M exp 


-^Phenom 




(MeV) 


(MeV) 


(e 2 fm 4 ) 


(e 2 fm 4 ) 


(fm) 








3U S 


2.79 


2.21 


327 


320±40 


1.19 


1.05 


1.05 


0.93 ±0.20 


32 S 


2.94 


2.23 


294 


300±13 


1.14 


0.96 


0.96 


0.95 ±0.11 


34g 


2.65 


2.13 


256 


212±12 


0.95 


0.94 


1.04 


0.91 ±0.11 


36 S 


3.46 


3.29 


241 


96±26 


0.63 


0.64 


0.90 


1.13 ±0.27 


38 S 


2.19 


1.30 


325 


235±30 


0.98 


0.98 


1.44 


1.50 ±0.30 


40 g 


1.54 


0.89 


431 


334±36 


1.17 


1.03 


1.17 


1.25 ±0.25 



It is evident from this definition (see also p. 65 7 of Ref. [^9|) that these are the appropri- 
ate quantities to be folded with the NN interaction to obtain the inelastic proton-nucleus 
form factors as described in Sec. ^|. The explicit radial form of the QRPA transition 
density, for the transition to an excited state \k > of the multipolarity A, is 

A/£(r) = p\{r) = J2 Va(r)^(r) < P\\Y x \\a > {X« - Y^}{u a v p + (-l)VM .(43) 

a>fi 

The HF±BCS ground state fl4*l"|) and QRPA transition densities (fUD of the considered 
Sulfur isotopes are those obtained recently by the Orsay group |^T[ (see, e.g., the radial 



shape of these densities in Figs. 10 and 11 of Ref. pl|). In particular, the QRPA charge 
ground-state and transition densities for 32 > 34 S have been shown to reproduce reasonably 
well the experimental charge densities of these isotopes (measured directly by electron 
scattering). It should be noted that these transition densities have been previously used to 
obtain the inelastic S±p form factors using the JLM effective interaction |21[| . Therefore, 



the comparison of the results obtained for the Sulfur isotopic chain will lend crucial 
information on the validity of reaction models (see further discussion in Sec. |j). 

Among the integral properties of the QRPA transition densities, the most important 
ones are the reduced transition probabilities B(EX |) and the ratio of the transition 
moments M™/M^, as determined by Eq. fl20|). In a simple collective model P5| , ^6 



the oscillation of the homogeneous neutron-proton fluid is so-called purely isoscalar if 
M"/M^ = N/Z. Consequently, a significant deviation of M%/M^ from N/Z would indi- 
cate the degree of what we call the isovector mixing in the considered quadrupole excita- 
tions. We will discuss this as well as other structure effects in more details in the next 
Section. 



4. RESULTS AND DISCUSSION 
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4.1. Elastic scattering 

For any nucleon-nucleus system, the accurate measurement of the elastic scattering 
provides essential information for the determination of the nucleon-nucleus optical poten- 
tial, which is further used to generate the distorted waves for the calculation of inelastic 
scattering amplitude in the DWBA formalism. The considered elastic 30— 40 S+p scattering 
data at different energies have been analyzed in the present work with the optical poten- 
tial obtained from Eqs. fl36|) and (p7|). The real optical potential was obtained using the 
CDM3Y6 interaction and the ground-state densities given by the HF+BCS calculation, 
while the WS imaginary potential was parameterized using the CH89 global systematics 
[39(1 . The best-fit optical model parameters are given in Table From the results plotted 



in upper parts of Figs. [l|-[| one can see that our folding potential reproduces the elastic 
data quite well. Although the CH89 parameters are mainly based on the nucleon-nucleus 
elastic data for medium-mass targets, they seem to be able to provide a reasonable esti- 
mate for the Sulfur isotopes as well. In the 30 S+p and 40 S+p cases, the CH89 parameters 
of the WS imaginary part of the optical potential are very appropriate and only the 
renormalization factor Nr of the real folded potential was slightly adjusted by the OM 
search to fit the elastic data. For the 32_38 S isotopes, the strengths Wv, Ws were also 
adjusted by the OM fit and the calculated elastic cross sections agree very well with the 
measurement. The obtained OM parameters (Table |3|) have been used further in all the 




Figure 1. Elastic and inelastic 30 ' 32 S+p scattering data at E/A = 53 MeV |21| in comparison 
with the DWBA cross sections given by the elastic and inelastic potentials folded with the 
HF+BCS ground-state and QRPA transition densities, respectively. The cross sections given by 
the isoscalar potentials alone are plotted as dotted curves. 



We note that the recently measured elastic 30 ' 32 S+p data at E/A = 53 MeV 
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are quite accurate and have more data points compared to other cases. To show the 
reliability of the semi-microscopic optical potential obtained in the present approach, we 
have further plotted in Fig. |2| the elastic 30,32 S+p data and the calculated cross section in 
ratio to Rutherford cross section (because the actual quality of the OM fit then becomes 
more transparent). One can see that the oscillating structure of the elastic cross section 
becomes more pronounced if plotted as ratio to Rutherford, and our OM results give 
indeed a ver 
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Figure 2. Elastic 30 ' 32 S+p scattering data at E/A = 53 MeV and results of the present 
OM calculation plotted in ratio to corresponding Rutherford cross section. 

In recent years, various experimental studies of reaction cross sections have been done to 
probe the nuclear sizes, especially the existence of halo or neutron skin in unstable nuclei 
[ f47f . It is, therefore, of interest to present the calculated total reaction cross sections 
<tr for the considered cases (see Table §). From the relative difference Actr between the 
total reaction cross section given by our folding analysis and that given by the mass- and 
energy dependent systematics for stable isotopes (see Eq. (1) in Ref. Pffi), one finds that 
our results agree with the systematics within about 10%. Actr is only 5% for 32 S (a stable, 
N = Z isotope), and we do not see any enhancement in <7r for the proton-rich 30 S isotope 
that might indicate a proton halo in this nucleus. To provide crucial data on this point, 
measurements at large center-of-mass angles are called for. Note that the total reaction 
cross section is mainly determined by the imaginary part of the optical potential, and the 
large angle data, if available, would reduce the uncertainty in the absorption strength and 
a R values can be extracted more precisely. Among the Sulfur isotopes under study Actr 
is largest for 40 S (about 18%) which is clearly due to the neutron skin in this neutron-rich 
(unstable) isotope. The neutron (proton) r.m.s. radii calculated within HF+BCS for this 



S+p, E, ab =53A MeV 




• 32S 

o 30S 
— Folding 
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Table 3 

OM parameters [see Eqs. ( |36|) and (p7|)1 used in the folding analysis of the elastic proton scat- 
tering on Sulfur isotopes. The real optical potentials were folded with the CDM3Y6 interaction 
and ground-state densities given by the HF+BCS calculation. Starting parameters of the WS 
imaginary optical potential and Coulomb radius were taken from the global systematics CH89 
1 39]. Parameters which were fixed during the OM search are given in boldface. 



Target 


E/A 


Nr 


w v 


W s 


r w 


a w 


r c 








(MeV) 




(MeV) 


(MeV) 


(fm) 


(fm) 


(fm) 


(mb) 


(%) 


3u S 


53 


0.940 


5.347 


3.723 


1.195 


0.690 


1.275 


669.4 


12 


32 S 


53 


0.901 


4.037 


3.480 


1.198 


0.690 


1.275 


643.3 


5 


34 s 


30 


0.930 


3.949 


6.315 


1.200 


0.690 


1.274 


937.5 


12 


36 s 


28 


0.947 


1.074 


6.968 


1.203 


0.690 


1.273 


945.7 


11 


38 s 


39 


0.935 


4.315 


4.809 


1.205 


0.690 


1.273 


894.5 


10 


40 g 


30 


0.990 


2.571 


7.926 


1.207 


0.690 


1.272 


1079 
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Act? is the relative difference between the total reaction cross section given by our foldin 
analysis and that given by the mass- and energy dependent systematics for stable isotopes 



nucleus are 3.44 fm (3.24 fm), as reported in Table 2 of Ref. |H . 



In an OM study of the elastic scattering induced by a neutron-rich nucleus, the isospin 
dependence of the optical potential should become more significant. While this effect is 
negligible in 30 ' 32 S+p cases, it becomes more and more sizable in 36 > 38 > 40 s_|-p cases. From 
the prediction of the HF+BCS calculation for the neutron rich 36 > 38 > 40 S isotopes (see, e.g., 
Fig. 10 in Ref. 0]), the extra neutrons do not form the halo-like structure and they 
distribute evenly both in the interior and at the surface. As a result, the ground state 
neutron density is larger than the proton density at all radii. This difference determines 
the isovector part of the optical potential Vq V (E, R), i.e., the microscopic estimation of the 



Lane potential jijj. Our results show that the maximal difference in the elastic scattering 
cross section, caused by the isovector potential Vq V (E, R), is ranging from about 15% in 
36 S+p case up to about 30% in 40 S+p case (see upper parts of Figs. |3] and |j). Although 
the difference found is still within the experimental errors of the measured elastic cross 
section, this result indicates that the isovector component of the proton-nucleus optical 
potential might be probed in accurate measurements of the elastic scattering induced 
by neutron rich beams. Such experiments could be an alternative study to (p,n) charge 
exchange reaction (where the isovector component of the nucleon-nucleus optical potential 
can be separately tested |5Tj|.[5l|). 



4.2. Inelastic scattering 

The present folding analysis of inelastic scattering data has been performed within the 
standard DWBA using the transition form factors defined in Eqs. (j38|) and (|39|). Since 



the effective NN interaction has been fine-tuned in the HF studies of symmetric |T5 



and asymmetric [18fl nuclear matter, it is expected to be a quite reliable input for the 
folding calculation and the test of the QRPA transition densities in the folding + DWBA 
analysis. One can see in Table |2| that the proton parts of the QRPA transition densities 
for the lowest 2 + states in 30 > 32 S give B(E2) values very close to the experimental ones. 
The corresponding neutron transition densities were found to have about the same radial 
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Figure 3. The same as Fig. [j] but for elastic and inelastic 34 S+p and 36 S+p data at E/A = 30 
[|| and 28 MeV @, respectively. 



shape (see Fig. 11 in Ref. plf ) and they give, therefore, M^+/M|+ « N/Z. From the 
DWBA results plotted in Fig. g one can see that the inelastic form factor calculated from 
the QRPA transition densities describes very well the inelastic 30 ' 32 S+p scattering data 
measured recently at = 53 A MeV [21]. This confirms the microscopic structure of 



the lowest 2 + states in 30 > 32 S predicted by the QRPA. Since the difference between the 
proton and neutron transition densities is very small for these two nuclei (especially for 
32 S) the isovector component of the folded form factor is negligible and the DWBA cross 
section is determined mostly by the isoscalar form factor. In other words, the lowest 
quadrupole excitation in 30 ' 32 S is predominantly isoscalar, with the best-fit moment ratio 
M = (M£+/M% + )/(N/Z) « 1. 

We note that, given the accurate QRPA transition density (which has a proton part 
very close to the experimental charge transition density ||21|| ) and well tested CDM3Y6 
interaction [15), the new proton scattering data measured for the (N = Z) stable 32 S 
isotope turn out to be a very good test case for our folding approach. The fact that, 
without any scaling of the folded form factor, the DWBA cross section agrees perfectly 
with the inelastic 32 S+p data confirms the reliability of the present folding + DWBA 
approach in probing the structure of the lowest 2 + states in Sulfur isotopes. This result 
also shows that the inelastic form factor folded with the CDM3Y6 interaction (which 
has been fine-tuned in the HF studies) can enter the DWBA calculation without a re- 
normalization factor, i.e., Nr = 1 for the inelastic folded potential. As a result, the 
higher-order contribution from the DPP to the real inelastic form factor seems to be 
weaker than that to the real optical potential. 

The DWBA folding results are compared with the 34,36 S+p scattering data in Fig. |3|. 
With the predicted B(E2) value about 20% larger than the experimental data for 34 S (see 
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Table 0), the agreement between the DWBA cross section with the inelastic data is still 
reasonably good for the 34 S+p system. In 36 S, the B(E2) value given by the QRPA is 
roughly 2.5 times larger than the experimental one, and the DWBA cross section for the 
2 + excitation in 36 S consequently overestimates the data over the whole measured angular 
region, especially at large scattering angles. The damping of the lowest 2 + excitation in 
36 S, leading to a transition strength B(E2) cxp about two times smaller than that observed 
for the lowest 2 + states in neighboring Sulfur isotopes, cannot be explained within the 
current QRPA model. The neutron transition strength predicted by QRPA turns out to be 
quite weak and gives a moment ratio M = (MJ+ / 'M%+) / '(N j 'Z) 0.64. Several structure 
models together with the DWBA analysis of the same data using the form factor folded 
with the JLM effective interaction have also failed to describe this nucleus. Thus, 

36 S nucleus still remains a puzzle to theoreticians, and the structure models should start 
being able to reproduce its B(E2) value. We will see below that some disagreement 
with the measured angular distribution remains after renormalizing the QRPA (proton) 
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Figure 4. The same as Fig. [T] but for elastic and inelastic 38 S+p and 40 S+p data at E/A = 39 
||§ and 30 MeV |2jJ, respectively. 



The DWBA results for 38,40 S+p systems are compared with the scattering data in Fig. [|. 
Even though the predicted B(E2) value is about 40% larger than the experimental one 
in 38 S (see Table |2|), the DWBA cross section underestimates the inelastic 38 S+p data 
by about 30%. This indicates that the QRPA neutron transition strength is somewhat 
weaker than that required by the inelastic data. The predicted moment ratio M% + /M%+ 
for 38 S is close to N/Z (the same as in 32 S) and it could also indicate a lack of neutron 
transition strength in the QRPA results. For 40 S, the predicted B(E2) is about 30% larger 
than the experimental value and M^+/M^ + ~ N/Z (see Table |2]). However, it is difficult 
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to discuss about the proton and neutron transition strengths of the 2 + state in 40 S based 
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Figure 5. The same as Fig. || but for elastic and inelastic 36 S+p and 38 S+p data at E/A = 28 
[ p3[ and 39 MeV [^4|], respectively. Inelastic (transition) form factors folded with renormalized 
QRPA transition densities were used in the DWBA calculation (see details in Text). 



Clearly, one can discuss quantitatively the transition strengths only when the theoretical 
model gives good description of both the B(E2) cxp value and the measured inelastic 
angular distribution. For that purpose, we have further made the same folding + DWBA 
analysis of 34_40 S+p systems using an inelastic form factor obtained from renormalized 
QRPA transition densities, as done earlier by Khan et al. [21|. Namely, the proton part 



of the QRPA transition density is scaled to give the experimental B(E2) exp value and 
the neutron part is adjusted by the best DWBA fit to the inelastic scattering data. The 
moment ratio obtained in this way is denoted as M exp in Table |2|. Note that for 30 ' 32 S 
isotopes, the QRPA transition densities give good description to both the B(E2) cxp value 
and the measured inelastic scattering data, so that M exp = Mqrpa- Among other cases, 
the improved agreement between the DWBA cross section and the data is clearly seen 
in Fig. [5] for the 36,38 S+p systems. One can notice that, even after the QRPA transition 
density is renormalized, there remains some disagreement with the inelastic 36 S+p data 
at medium angles. The scattering at medium and large angles is known to probe the form 
factor at sub-surface distance, and this effect might indicate a deficiency in the radial shape 
of the QRPA transition density for 36 S which cannot be changed by the renormalization 
procedure. The damping of the 2 + excitation in 36 S can be judged by the scaling factors 
of the QRPA transition density: after the proton transition density is scaled down by a 
factor of 0.63 to give B(E2) = B(E2) exp « 96 e 2 fm 4 , the neutron transition density needs 
to be renormalized by a factor of 0.88 to fit the inelastic data. Although the obtained 
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moment ratio (M exp ~ 0.90) is larger than that predicted by QRPA for 36 S, it still shows 
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Figure 6. 38 S+p (real) optical potential (left panel) and inelastic form factor for the 2 + exci- 
tation (right panel) folded with the HF+BCS ground-state and renormalized QRPA transition 
densities. Isoscalar and isovector components of the potentials are plotted as dashed and dotted 
curves, respectively. 



The results obtained for the 38 S+p system are quite interesting. After the proton 
transition density was scaled by a factor of 0.85 to give the experimental B(E2) exp value, 
the neutron transition density needed to be renormalized by a factor of 1.25 for the best 
agreement of DWBA cross section with the inelastic data. As a result, the difference 
between the proton and neutron transition densities becomes larger and enhances the 
isovector contribution to the DWBA cross section. The 2 + cross section given by the 
isovector part of the 38 S+p form factor is about 25-30% of the total 2 + cross section. Such 
a significant contribution is due to a strong neutron transition strength (M£+) compared 
with the proton one (M^ + ) which gives the ratio M exp = 1.44. The isovector mixing is 
slightly weaker in the elastic channel, with the cross section given by the isovector part of 
the 38 S+p optical potential of about 20% of the total elastic cross section. In Fig. one 
can see that the difference (p p — p n ) between the proton and neutron transition densities 
leads to an enhancement of the isovector part of the folded 2 + form factor [see also 
Eqs. ( p6[ ) and (|3*0|) 1 at distances around R = 4 fm, where amounts up to 15% of the 
total form factor strength. Since the inelastic cross section is directly proportional to the 
square of the form factor, the 2 + cross section given by the isovector part of the folded 
form factor alone (see Fig. |5[) becomes around 25-30% of the total 2 + cross section. For 
the elastic scattering the difference (p p — p n ) between the proton and neutron ground- 
state densities is smaller, and the isovector part of the folded optical potential or the Lane 



21 



potential is about 7-9% of the total potential strength at radii near R = 4 fm (see left 
panel of Fig. which leads to an isovector mixing of about 20% in the elastic 38 S+p 
cross section. 

For 40 S+p system, the isovector mixing becomes stronger in the elastic channel, with 
the cross section given by the isovector part of the 40 S+p folded optical potential of about 
30% of the total elastic cross section. However, the situation in the 2 + inelastic channel 
is still uncertain because the best-fit ratio M exp = 1.17 is based on two data points only. 
One clearly needs more experimental information in order to have a reliable estimate for 
the transition moment ratio in S. 

The comparison of the results using the same structure input but with two different 
reaction models (JLM approach and the present folding model) is of special interest, 
because of the extensive (p,p') data available along the Sulfur isotopic chain, including 
elastic and inelastic angular distributions. Although the JLM approach gives, in general, 
good agreement with the data, it appears to have some difficulties to reproduce both 
the first minimum and the elastic distribution at large angles (see Figs. 13, 14, 15 of 
Ref. pifl , Fig. 9 of Ref. [^5|] and Fig. 7 of Ref. []52|). It is interesting to note that this 



feature does not happen with the present folding calculations. For the considered 2 + 
excitations, the DWBA cross sections given by the inelastic form factors folded with the 
JLM interaction seem to overestimate the inelastic data in most cases which results in 
the smaller transition moment ratios compared to our study. First, one should note the 
radically different approaches of these two models. In the JLM approach, one uses the local 
density approximation (LDA) to obtain the local optical potential from that calculated for 
infinite nuclear matter, while in the folding model both the optical potential and inelastic 
form factor are calculated in the Hartree-Fock manner (with antisymmetrization treated 
properly) using a fine-tuned density dependent CDM3Y6 interaction that is based on the 
G-matrix elements of the Paris NN potential for finite nuclei. Therefore, the mentioned 
problem might well be due to the validity of the LDA and suggests further investigations to 
improve the JLM parameters. Another clue is that a better description of elastic data by 
the folding potential might also be due to the fact that the imaginary optical potential is 
based on CH89 global systematics fl59| , with its strength slightly adjusted to fit the elastic 
data. In this sense, the imaginary WS (elastic and transition) potentials obtained in our 
approach could serve as a good reference to improve the microscopic description of the 
imaginary part of the optical potential and the transition form factor for the considered 
cases. As discussed above, the 36 S case remains a common problem for both reaction 
models and the deficiency should be due to the structure model. Finally we note that the 
best-fit transition moment ratios deduced from our folding + DWBA analysis turned out 
to be quite close to those deduced earlier (see M exp and Mphenom values in Table ||) from 
the DWBA analysis of essentially the same data |T],|5| using standard collective model 
form factor. This is very satisfactory in the sense that the microscopic analysis of the 
properties of stable and exotic nuclei, which is the final goal of this kind of approach, is 
independent from the reaction model, provided this reaction model is accurate enough. 
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5. CONCLUSION 

A consistent single folding formalism for the proton-nucleus optical potential and in- 
elastic form factor, using a realistic density dependent CDM3Y6 interaction based on 
the G-matrix elements of the Paris NN potential and the microscopic ground-state and 
transition densities given by the HF+BCS and QRPA calculations, respectively, has been 
applied to study elastic and inelastic scattering of 30_40 S isotopes on proton target. The 
contribution of the isovector part of the folded optical potential and inelastic form factor 
to the calculated elastic and inelastic cross section has been considered explicitly in all 
cases. 

For the proton rich 32 S nucleus, the results of the OM and DWBA analyses reproduce 
very well the measured elastic and inelastic scattering data and do not indicate the ex- 
istence of a proton halo in this nucleus. Our results have shown quite a strong isovector 
mixing in the elastic and inelastic 38 ' 40 S+p scattering. In particular, the best fit strength 
of the isovector part of the form factor in the inelastic 2 + channel for 38 S is significantly 
stronger than that predicted by the QRPA. The detailed folding + DWBA analysis of the 
considered data has closely reproduced the systematics found earlier (based on standard 
collective model form factors) |pljp5| for the neutron-proton ratio of the 2 + transition 



moments in neutron rich Sulfur isotopes. 

In conclusion, the present folding + DWBA approach has been proven to be accurate 
and reliable in extracting the neutron and proton transition strengths M™^ of the nuclear 
excitation induced by the inelastic proton scattering. Moreover, our approach is very 
efficient in testing the microscopic nuclear densities and determining the isospin character 
of the considered excitation (in structure of the excited state as well as in the inelastic 
angular distribution). 

These results confirm again that proton scattering induced by the neutron rich beams 
remains a very effective tool in nuclear structure study. Such measurements would allow us 
not only to test the ingredients of the microscopic structure models (like QRPA) but also 
the isospin dependence of the proton-nucleus optical potential and inelastic form factor 
(prototypes of the Lane potential) and, consequently, the isospin dependence of the in- 
medium NN interaction (which is a key to specify the equation of state for asymmetric 
nuclear matter fff31). 
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